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I. INTRODUCTION 



Until the pioneering work of Drukarev and Levin on the use of QCD sum rules in systems of finite nucleon density 
, descriptions of nucleons in nuclear matter had mostly been based on models of N-N interactions or the direct 
use of meson exchange potentials in quantum hadrodynamics ||. These models have serious problems at both the 
short-distance and long-distance scales. At short distance one must treat the complex structure of hadrons, for which 
we now believe that QCD gives the correct description. At long-distance scales it seems that meson cloud effects 
must be included, and neither quantum hadrodynamics nor current treatments of QCD are adequate; while chiral 
perturbation theory is very successful in a quantitative treatment of pionic effects at low energies. Recently a 
model was developed for use in QCD sum rules with a nucleon current that has a Goldstone boson component, which 
has the main features of both microscopic QCD and chiral perturbation theory Q . In the present paper we extend 
this model for the study of nucleons in finite nuclei. 

The strongest experimental evidence for the utility of introducing an explicit meson cloud component for the nucleon 
is the data on sea quark distributions obtained from Drell-Yan experiments. The NMC/CERN experiments |?J have 
given evidence for violation of the Gottfried sum rule 

f — [F* (x) - F? (*)] = \f { [ul (x) - d; (a)] + 2 [u* (x) - d; (x)} } dx. (1.1) 
Jo x 6 Jo 

If the up-type and down-type sea quark distributions are equal, as would be expected from QCD, then the value of 
this integral is 1/3. The NMC result is 0.24 ±0.016 at Q 2 = 5 GeV 2 over the interval 0.004 < x < 0.8. Similar effects 
were found in Drell-Yan measurements from Fermilab/E866 jq|, where the ratio of down to up sea quark distributions 
is much larger than 1.0. Both of these results differ from the predictions of perturbative || and nonperturbativc |nj 
QCD calculations. One approach to account for this discrepancy in sea quark distributions in hadronic models is to 
include a meson cloud for the nucleon |ll] [l| . Our approach is to treat the pion as a basic Goldstone boson field at 
low energy, along with quarks and gluons, with separation of long-distance from medium- and short-distance effects. 

Another strong motivation for introducing a correlator with a Goldstone boson component for use in QCD sum 
rules is given by chiral perturbaton theory. It was shown almost two decades ago [ ff4[ that the leading nonanalytic 
contribution (LNAC) to the nucleon mass in chiral perturbation theory is a m^ correction, with the lowest-order chiral 
logs not contributing. In the standard formulation of QCD sum rules one does not obtain such a m^ term jL5|. In 
Rcf H it was shown not only that one obtains the m^ LNAC, but that the known value of the mass shift (15 MeV) 
can be used to help determine the parameters of the model, in spite of the fact that the sum rule method is only 
accurate to about 10% for the nucleon mass. It is not possible to obtain the LNAC m^ term with the usual quark 
fields correlator We make use of this result in the present work. 

It should also be noted that in the QCD sum rule method the microscopic QCD calculation gives m 2 ln(m 2 ) terms 
for the nucleon in vacuum and also pm, terms for nucleons in the nuclear medium, both in contradiction to chiral 
perturbation theory. It was shown that both for the vacuum | ft5l and for nuclear matter |]l6| a careful treatment of 
the contunuum part of the phenomenological correlator enables one to cancel these terms, which are inconsistent with 
chiral perturbation theory. We made use of this cancellation of unwanted logs in the formulation of our model and 
use this mechanism in the present work. 

We review Goldstone boson/QCD model of Ref. Q for a nucleon in the vacuum in Section 2. In Section 3 the 
surface effects on the pion propagator are quantified using a p-wave ttA optical potential where the A(1232) resonance 
is dominant. We extend the model for the microscopic correlator from zero to finite density in Section 4, with the 
goal of calculating nucleon mass shifts due to the pion cloud effects in small symmetric finite nuclei. The sum rules 
are then constructed and evaluated in the final sections. 



II. REVIEW OF NUCLEON CORRELATOR WITH PION CLOUD IN VACUUM 

In this section we review the model for the nucleon correlator with explicit Goldstone boson fields developed in Ref 
. The nucleon current has the form 

V N (x) = c lV N ^{x) + c 2 V N ' % (x) (2.1) 

where the constants Cj are the amplitudes of the composite field operator without and with the Goldstone boson field. 
This operator is used to construct the nucleon correlator. We ignore possible contributions to the nucleon from the 
strange part of the meson cloud. The proton field operator without the pion cloud is chosen to be the current jl7) 

rf>°( x ) = e abc [u a (x) T Cjy(x)] ^-fd^x), (2.2) 



2 



where a,b,c label the color indices, and u(x),d(x) are the up and down quark fields. C is the charge conjugation 
operator. The lowest-energy contribution to the phenomenological dispersion relation for the correlator, the pole term 
giving the proton intermediate state, depends on the transition matrix element 

(0|rf' (p)|proton) = X p v(p), (2-3) 

where X p is a structure constant and the Dirac spinor v(p) is normalized by 

v(p)v(p) = 2m. (2.4) 



The neutron current r] n '°(x) can be obtained by the interchange of up and down quark fields in Equation (2.2). The 
nucleon current including the explicit Goldstone boson field is taken to be 

rf>*(x) = ^dMx) ■ Tf-f^x), (2.5) 

K 

where <fr„(x) is a massless pion field, r is the I-spin operator and A^ is a D = 1 scale parameter. The structure 
constant of this current is defined by the relation 

(0 \rf'*(x) | proton) = X' p v(p). (2.6) 

The p-wave coupling of the pion to the nucleon current and experience with hybrid mesons jOI suggests that X' p 2 <C Ap, 
so we expect little contribution to the phenomenological expression for our nucleon correlator from the resulting pole 
due to its very small residue. 

The full correlator for the proton with the pion cloud is 

W(p) = i [ d A xe lpx (0 \T [rf(x)ff(0)]\ 0) 



d 4 xe tp - x (0 \T [rf'°(x)ff'°(0)] | 0) + (1 - c?)i J d A xe lp ' x (0 \T [rf' w (x)ff ^(0)]] 0) 

= c\Xl^ \p) + (1 - cl)n^\p). (2.7) 

With the appropriate changes of current and isospin considerations, the neutron correlator is obtained. This model 
focuses on the long-range effects of the pion cloud and therefore docs not include pion-quark interactions and the 
coupling between currents that include the pion cloud and currents that do not. It is assumed that these short-range 
contributions are already accounted for in QCD through the condensates, and their inclusion in the model would 
result in double counting some of the pionic effects on the nucleon correlator. 

The part of the proton correlator without the pion cloud can be found in Ref. |fj"7|| . The nucleon correlator with the 
meson cloud for the microscopic evaluation of the sum rule has the form 

U M(x) = Dl a0 (xh a n ip -°Hxh , (2.8) 

where n^ p,0 ^(x) is the coordinate space nucleon correlator without the pion cloud |l7j ] and in the chiral limit DJ Q(3 (x) 
has the form 

D i^ x) = - • (2 - 9) 

In momentum space the nucleon correlator with the pion cloud is 

n^)(p) = -2 j^±zM^mtJi (2 . 10a) 



(27T)* (p-kY 

= Mdd ) b) + n^ n ) ( P ). (2.10b) 
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(a) (b) (c) 




FIG. 1. The diagrams corresponding to the pion cloud contributions to the nucleon correlator in vacuum. The solid lines 
correspond to quarks, the helical lines to gluons and the dashed lines to Goldstone bosons. The circles and ovals correspond to 
the vacuum condensates. 

The pion cloud contributions to the nucleon correlator are shown in the diagrams of Figure [l], and in momentum 
space are 

ni fc (;p) 
n w (p) 
ni e Op) 

Details are given in Ref. [§. The usual QCD sum rule methods are used, as in Ref. [jITj, but with a modification of 
the treatment of the continuum given by the work of Ref. |isj| . In addition to the calculation of the nucleon mass 
the magnetic dipole moments of the nucleons were calculated. From the comparison of the results for the magnetic 
dipole moment to the result without the pion cloud component flil it was estimated that the pion cloud component 
has roughly the same probability as the cloudless component, or c\ ~ c\ ~ 0.5. 



^-^/M-P 2 )t (2.Ha) 



2 

P 2 \n{-p z )i> (2.11b) 



i (gg) 2 , . i 
2 3 3^ 2 Aj 

T -/ln(V)^ (2.11c) 



2 n 3^ 6 Aj 

■ 2 rr . nn\ 

\n{-p 2 )j>. (2. lid) 



-i (qq) (q9s<7 ■ Gq) , , .. , 



2 5 ^ 2 Ai 



FIG. 2. The diagrams corresponding to the lowest-dimension chiral corrections to the nucleon correlator in vacuum 

An important observation is that the diagrams shown in Figure ^ gives the main LNAC m\ terms. The largest 
contribution is from the first term shown in Figure ^a, which was shown to give the following contribution to the 
nucleon mass 



AM p 



em%g\ 



160/3 2 /2 

-20 MeV, (2.12) 



with the choice of A2 = Mp/g^jy, which is quite reasonable. 

The contribution to the nucleon mass from the pion cloud was found to be negligible, consistent with the assumption 
X'p <C Ap. We use this value for A^ in the present work. This completes the model within the uncertainty of the 
probability of the meson cloud component. 
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III. PION PROPAGATOR IN FINITE NUCLEI 



Extending the model for the nucleon correlator with the pion cloud from vacuum to finite nuclei requires an 
appropriate modification of the pseudo-Goldstone boson field in medium. In order to be able to able to apply the 
model at high momentum transfer, which is needed for various applications, we make use of conventional multiple 
scattering formalism [see, e.g., Ref. pc| for a review] rather than chiral perturbation theory. This is a mean field 
method in which the pion propagator in the nucleon medium is given by the pion optical potential. 

The propagator for a free massless pion has the form 

m) = (3-D 

In the medium the pion acquires an effective mass given by the optical potential. The propagator for a pion with 
energy fco i n & nucleus of mass number A is 

D\(k) = \ (3.2) 

k 2 + Jl(k ,k;p,A) + ie 

where H(k , k; p, A) is the pion self-energy in the medium and p is the baryon density of the matter. The pion self 
energy is related to the pion- nucleus optical potential V op t(ko) by 

n(fc , k, k') = -2k (k' |C P f (to) | *) , (3.3) 

where k and k' are respectively the incoming and outgoing three-momenta of the pion. We begin constructing the 
n A optical potential by looking at the partial wave expansion of the irN scattering amplitude 



T(k',k)=J2Pi fe [(i + 1) /7,i+H+i/i,i-(w)]fl(oosfl)- iff ^ [// 

i V i i 



l+(«)-/i,l-(w)]iY(c 



(3.4) 



The partial wave amplitude is related to the S-matrix by 



fa = -4- (SaM - 1) (3-5) 

2i|fc| 

where the scattering matrix is found using the experimentally obtained phase shifts S a (u>) as 

S a (w) = e 2i5 "K (3.6) 

In this notation a indexes the isospin /, orbital angular momentum I and total angular momentum J of the partial 
wave. The phase shifts are expected to depend on momentum proportional to k 2l+1 , and we define the constants a a 
as 

Sa (3.7) 



k 2l + l 



These constants are the so-called scattering lengths (027) for s- waves and the scattering volumes (62/, 2 j) for p- waves. 
These parameters for d-wave and higher order partial waves are very small and have negligible effect on irN physics, 
so we truncate the expansion of the scattering amplitude at I = 1. Then the irN scattering amplitude can take the 
form 

F{h', k) - a(fc ) + ar(k ) (f • f) + [b(k ) + &j(fc ) {t- ¥)] k' ■ k + i [b sf (k ) + b Isf {k ) (t ■ f)] S ■ (%' x fc) , (3.8) 

where t and f are the isospin matrices for the pion and the nucleon. The energy-dependent coefficients can be related 
to the scattering volumes and lengths at threshold; for example, b = -| (4633 + 2613 + 2631 + b\\). Because we are 
concerned in this work with spin/isospin symmetric nuclei, the isospin dependent and spin-flip parts of the scattering 
amplitude vanish. Furthermore, the p-wave contibution is considerably larger than the s-wave contribution at low 
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and medium energies because of the dominance of the A(1232) 3, 3 resonance in ttN scattering. We are then left with 
the simple scattering amplitude for the spin and isospin averaged scattering of a pion off of a nucleon, 



F(k',k) = b(k )k' -k. (3.9) 

In order to take into consideration the off-energy-shell effects in our optical potential, we modify the result of ( 
by considering the p-wave t-matrix for the ttN scattering in the separable form 

(k' |^~ wavc (k )\ k) = -b{k )g{k')g{k)k' ■ k, (3.10) 

where g(k) is a form-factor that roughly parametrizes the inelastic channels in the nN scattering at high energies. 
g(k) is chosen be unity when k is at the on-energy-shell value, and to approach zero as k gets large, where k = 
We use a monopolc paramctrization of the form factor 

\k\ 2 + A 2 

with vertex factor A taken to be 700 MeV, for the description of the off-shell properties of the amplitude. 

The energy dependence of the function b(ko) is modelled by using the Breit-Wigner form for the A(1232), giving 

6(fco, k) = *°1 ( - h , v (3.12) 

V 3?rm^ J iv A - k - § r A (k, ko) 

where ua is the energy of the 3, 3 resonance and / is the pseudovector ttN coupling constant. We model the resonance 
so that the three-momentum dependence of the width Ta is cut off at low momenta and at the momentum where the 
resonance is saturated: 

Ta (fc, ko) = { k 3 k m <k<k A . (3.13) 

v ' 3 47r m£ ko 



The lower limit starting at k m = 30 MeV ensures that the A(1232) does not vanish even at very low pion momenta, 
and the upper limit at the three-momentum k A where the Breit-Wigner curve peaks provides a mechanism for the 
resonance to disappear at high energies as expected. 

To obtain the pion- nucleus optical potential for elastic scattering we make use of three approximations: 

• low density 

• impulse 

• local density. 

The low density and impulse approximations lead to the standard first-order in p form for the optical potential, and 
the local density approximation is used for finite nuclei to treat the interaction at the nuclear surface as a function of 
p. The p-wave pion optical potential for elastic scattering by an A-nuclcon nucleus can then be written to first order 
in density as 




>' |C p f| k)=A (k' ICIH k) P [k k') . (3.14) 
p{k — k') is the Fourier transform of the ground state nuclear density p(r). The resulting optical potential in coordinate 



space, using Eq.(3.1C), has the form 

V(f) = Ac(E^) V -P(r) V ■ (3-15) 

This is the so-called Kisslinger potential [^l| without the s-wave scattering term. It has been shown to give an accurate 
treatment of medium-energy pion- nucleus elastic scattring p0[ . Due to the p-wave nature of the interaction derived 
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from the A(1232) resonance the interaction takes place mainly on the surface, and indeed vanishes in the interior of 
the nucleus if a constant density is used. This will be seen to be an important aspect of our present work. 
We model the spatial dependence of the nuclear density as a second order harmonic oscillator 



p[r) = po 1 + a - 



exp 



(3.16) 



where po is the density of nuclear matter, c is a size parameter, and a provides some control over the shape, allowing 
for modelling of heavier nuclei (like 40 Ca) and more complicated distributions (for example, 16 with a local minimum 
at the center). 

We now turn our attention to the pion propagator in medium. The propagator (|3.2|) is expanded in density as 



D\{k) = [D\{k)] po=0 + Y,~P^ 



S 



Sp(q) 



po=0 



2! 



91 12 



Sp(qi) Sp(q 2 ) 



(3.17) 



Po=0 



where the first term is just the free pion propagator (3.1) and q is the momentum transfer in the scattering. The 
scattering angles are averaged over, and the pion propagator in the nucleus takes the form 



Dm = Dj(k) + 



^ 2 {A-l) Pr> f 2 uo A 



g (k) Fi (fc, a, c) [uj a - k - -r A (k, k ] 



om 2 



O (pg) , (3.18) 



with the nuclear structure information contained in the function 



F t (k.n.r) = [ 1 + y 



c 2 k 2 



2c 2 k 2 



c 2 k 2 



c 2 k 2 



2c 4 k 4 + Ac z k 



.212 



c 2 k^ 



exp (-c 2 fc 2 ) 
exp (-c 2 P 



(3.19) 



IV. THE NUCLEON CORRELATOR WITH PION CLOUD IN MEDIUM 



To construct the in-medium nucleon correlator with the meson field, we start by rewriting the vacuum matrix 
element of the product rj(x)fj(x) in Eq.( |T7]) as a nuclear matrix element so that the proton current correlator takes 
the form 

D* (p) = i ( irtee* - " (A \T [rf (x)fj p (0)]| A) 



(4.1) 



= ^ (A |T | A) + cl.J Ae- (A \T [rf' w (x)ff' w (0)]\A) 



where \A) is the ground state of the finite nucleus. The part of the in-medium proton correlator with the pseudo- 
Goldstone field is 



S b u b ' A (xhvCS^ T (x)C^ 



+ - [d a d p D* A {x)) la rS~ A (x)l u l^ S bb A (x) lv CS^ r (x)Cy, 



[4-quark terms] 



1 



This is simply the expression for the meson cloud part of the correlator, given in Ref. j6|, where the quark propagator 
and pion propagator in vacuum, S^ b (x) ani 
in the nucleus, respectively. For example, 



r (n,0), 



«/3 



(4.2) 



and pion propagator in vacuum, S® (x) and DJ(x), are replaced by S^ b A (x) and D A (x), the quark and pion propagators 



S* b A (x) = (A\T[q a (x)q b (0)]\A) 



(4.3) 
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D\ a p(x) is the nuclear analogue to Eq.( |2.9| ), replacing the free pion propagator D*(x) with D\{x). As with the case 



in vacuum, the neutron correlator IL^'°'(x) is found simply by an interchange of up and down quark fields in the 
current rj(x). With a spin/isospin symmetric nucleus in the chiral limit, the in- medium correlator with a pion cloud 
has the form analogous to Eq.(|2.8|), 



For simplicity we work in momentum space, where the in-medium proton correlator has the form 



r(P.O), 



n>) = c ?n<f°>(p)-(i 



d 4 k 

(2^ 



Di(k)m { r ] (p-w 



where we made use of the Fourier transform of D\ a g{x) : 

D\ a0 (p) = -PaPffDHp). 



(4.4) 



(4.5) 



(4.6) 



Since the nucleus has a four-momentum u M , breaking Lorentz invariance, there are additional Dirac structures in the 
nucleon correlator. Under the assumption that the nuclear ground state is invariant under parity and time-reversal, 
it was shown [fij that the in-medium nucleon correlator is 



if; (p) = n ^ ( P 2 , P ■ u) + n£> { v \p- U )i> + n ™ ( P 2 , P -u 



r(«), 



(4.7) 



The individual scalar functions LT^ (p 2 ,p ■ u) depend only on the Lorentz invariants p 2 and p ■ u and can be projected 
out of the full nucleon correlator using the following formulae : 



U^(p 2 ,p-u) 
rtv\p 2 ,p.u) = 



<y,p- U ) 



i 



Tr[II»] 
1 

2 — (p ■ u) 2 
1 

p 2 — (p ■ u) 2 



iTr^(p)]-^Tr[^n»] 



(4.8) 



p ■ u 



Tr[^lF»] 



We define in-medium matrix elements of an operator O (x) as an expansion in local nuclear density p 

(O(x)) p =(A\O(x)\A)-(0\O(x)\0) 
= p(N\0(x)\N) + ---, 



(4.9) 



where (N \ (x) \ N) is the single- nucleon matrix element of O (x). The in- medium quark condensate is related to the 
ir-N sigma commutator by jjj 



VttN 



m d 



(4.10) 



We take a^N — 47 MeV p!| . The breaking of the Lorentz invariance of the medium by the introduction of the nuclear 
four-momentum leads to condensates with additional Dirac structures . The dimension three in-medium vector 
quark condensate, with the appropriate color and flavor weighting, is proportional to the nuclear density 



(4.11) 



The in-medium gluon condensate is found via the trace anomaly [E4 25 to be 

(g 2 G^G^) p = - p 2 6 39tt 2 MeV. 



(4.12) 



The values of the five-dimensional mixed quark-gluon condensates in medium are not well known, although significant 
cancellation between diagrams containing these matrix elements occurs much like similar terms in the free nucleon 
correlator in the chiral limit. Furthermore, it was shown in p2fl that the nucleon sum rules are insensitive to the 
values of these condensates, so we ignore terms proportional to the mixed condensates in our sum rule analysis. 
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The nucleon correlator is analyzed in the rest frame of the nucleus, where the nucleus four-momentum is u = 
(M, 0, 0, 0) with the mass of the nucleus M « Am (neglecting the average binding energy per nucleon in the nucleus 
as small compared to the nucleon mass m). In this frame p ■ u — 2AIpo, where po is the energy of the interpolating 
field. We carry out Borel transforms with respect to — p 2 , using the relationship 

(-J+p) =4m 2 . (4.13) 

In the nuclear rest frame, po — ^- (3m 2 — p 2 ), so when —p 2 — > oo, p ■ u — > oo while the ratio of the two invariants, 
analogous to the Bjorken scaling variable in DIS, remains fixed and the OPE is applicable. 

We restrict our calculation of the nucleon correlator to first order in the nuclear density. We can then identify three 
distinct types of quark current diagrams that contribute to the correlator: diagrams that contribute to the part of 
the in-medium correlator without the pseudo-Goldstone field (Figure ||) , diagrams that include the pion field where 
the quarks and gluons interact with the medium and the pion maintains its in- vacuum properties (Figure ^), and 
diagrams where the quarks and glue interact with the QCD vacuum and the pion scatters in the nucleus (Figure |^). 
The shaded boxes in these diagrams represent interaction of the field with the nuclear medium to first order in the 
density expansion. Diagrams where the quarks and pions both interact with the medium, or where the pion scattering 
includes two-body correlations would come in at second order in density and are not included in this work. 




FIG. 3. The diagrams corresponding to the contributions to the nucleon correlator in medium without the pion cloud. The 
shaded boxes represent the interaction of the fields with the medium to first order in density. 




FIG. 4. The diagrams corresponding to the free pion cloud contributions to the nucleon correlator in medium. 








FIG. 5. The diagrams corresponding to the in-medium pion cloud contributions to the nucleon correlator in medium. 



The first two sets of diagrams are most easily calculated in coordinate space using (4.4), after which they are Fourier 
transformed into momentum space. The contributions at first order in the density from diagrams without the pion 
cloud are found to be [fp^] 



n (2a,p)(p) 



n 



(2b.p) 



n (2c,p)(p) 

n 



n 



(2e+2/,p) 



(P) 



n 



(2i,p) 



(?) 



t^W<i) p [Hp 2 +i>u-p] in (-p 2 



(2ir 



2 7 7T 4 

3 

i 



(g 2 G-G)J\n(-p 2 ) 



<2,i, P )(p) = -^(m) {q f q) p ^f 



3 2 2tt 2 

i 



(W) (g 2 G-G)+(qq) p (g 2 G-G) 



u (^ P )(p) = Wq) p (g 2 G ■ g) 

2 2 i 1 

R(2h,p)(p) = — — (qq) {qq) P i>^ 



,1 ,u • p 
hi- - 2j>—f 



3 

2 2 m 



9 2 s (qq) (qq) 



(4.14) 



Here we have used a factorization approximation to evaluate the effects of the four-quark and six-quark condensates, 
where to first order in density we take (qTqqTq) p « (qq) (qq) p and (qTqqVqqYq) p w (qq) 2 (qq) p , respectively. While the 
factorization of these condensates is unjustified in medium [ p6| , our primary concern is a comparison of the nucleon 
mass shift due to the pion cloud terms and the shift due to terms without the pion cloud, so we find the approximation 
acceptable for the problem at hand. 

The diagrams in Figure ^, where there are first-order interactions of the quarks and gluons with the medium but 
the pions interact only with the vacuum, provide the following contributions at first order in the density: 



n (3a, P )(p) 

n (3c,p)(p) 



n 



(39, p) 



(p) 



-TTT Wq) n [ip 2 + ^i> u ■ p] p 4 ln (-p 2 ) 



2 6 3% 

n^(9 2 G-G) p M^(-P 
— 



H-(3h,p)(P) = 



2 8 3 3 tt 4 A 4 . 
i 



4T4 Wl) P (alG ■ G) [2 ■ hip 2 + lljm-p] In (~p 2 ) 



qq) (qq) p i>p W-p ) • 



2 2 3tt 2 A 4 

All diagrams which would be expected to contribute to scalar structure of the nucleon correlator vanish. 



(4.15) 
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The calculation of the terms represented in Figure || where the pions interact with the nucleus and the quarks 
and glue interact within the QCD vacuum are more difficult to calculate. Instead of calculating these diagrams in 
coordinate space as we did with the other in-mediu m di agrams, which would require first a rather daunting Fourier 



transform of the in-medium pion propagator (Eq. 3.18 ) and then an inverse transform of the resulting expression, 
we perform the calculations in momentum space using Equation ( |4.5| ) as the starting point. For example, the term 
corresponding to Figure gives the linear-density contribution 



(i f - - 2 
n (4 a,p)(p) = 2 i 07r 8 A 4 J dk {aw) F 1 (k,a,c) 

[„ , (loa - fa) fco + fr A (fc) (p - fc ) 4 ln [-(p- k Y 

X / dfc fc - " 



(wa - fco) 2 fc 2 + i (r A (fc)) 2 (fc 2 - fc 2 + ^) 2 



[(2p -k~k 2 )ft- k 2 p] 



where £ = 87r ^ ^ m Vg°'^ ,r " A ■ The Borel transform with respect to — p 2 , B, eliminates terms proportional to p ■ k in 

the expansion of (p — fc) 2 ™ ln[— (p — fc) 2 ], simplifying the integrand so that the term proportional to ft vanishes from 
symmetry in the angle-averaging. We are left with integrals over the pion energy fco and the magnitude of the pion 
three-momentum fc. The integrand also picks up a Borel dependence of the form Af|™ exp(— fc 2 /A7|). The energy 
integration can be completed analytically, where we find that the factor exp(— k 2 /M^) — » 1 and the dependence on 
the Borel mass can be taken out of the integral. This first diagram is found to be 



£n (4a , p) (p) = (a, c) -^fj>M% (4-16) 

A*. (271") 



where 7(a,c) is the remaining one-dimensional integral over the pion three-momentum. Defining T(fc) = fcor(fc,fco), 
the integral takes the form 

, , ~i f°° - r- w LO A k - k 2 + (k) 

1{a > c) = 2^ «( fc >«> c )(5( fc )) — 2 2 ,;i N 2 - ( 4 -17) 



( WA fc-fc2) 2 + l(f (fc))' 



We calculate this integral numerically. 

The other diagrams shown in Figure || yield the following contributions to the nucleon correlator: 



<i T , v 4tt 
B^{Ac,p)ij>) = tjI {a, c) 2 4 J {qqf 

7r d 



^n (4 b !p) (p) = -— 7 (a, c) — — F (qq) M B 



7}n (4e , p) (p) = -|i7 (a, c) -Lj* < ff 2 G • G) M 2 

(i 1 1 

£II(4/,p)0) = -^4 7 ( a > c ) ^4^7^ (W> (l9s°--Gq) -^gj 

(p) = -|-7 (a, c) (ft) (g 2 G ■ G) 



^7, „,,„(/>) = ^7(a,c) ^^(ft) 3 ^, ,-J.Im 



where 7 (a, c) in each of these expressions is the same as Equation (4.17) 



V. CONSTRUCTION OF THE IN-MEDIUM SUM RULES 
A. In the Chiral Limit 

The phenomenological ( "right-hand" ) side of the sum rules, is expressed as a dispersion relation with a lowest-lying 
hadronic pole of in-medium mass m* and residue A* 2 clearly separated from a continuum of higher-energy states with 
threshold s* §: 
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TJP fV> _ \*2„2 P 



- / ds 

7T 



(s,u- (?) 



P 



(5.1) 



Guided by the vacuum calculations, we assume that the structure parameter associated with the current including 
the pion cloud X* 2 <C A* 2 , so there is no additional pole considered in our model of the phenomenological correlator. 
Rather than find the values of the parameters in medium, we follow the appoach of Drukarcv and Levin and determine 
the shifts of the in-medium parameters from their vacuum values jjj. We remove the contributions to the nucleon 
correlator due to the pure in-vacua diagrams by defining the quantity 



(i) 



(5.2) 



where i labels the three structures in the sum rule (1, and -jl, respectively). On the phenomenological side of the 
sum rule we can express (j>^ in terms of the parameter shifts and the Borel transform of the RHS in vacuum as 



[M%, a) 



Am— + AX 2 — t— + As — 

dm dX 2 osq 



BU 



0,RHS 



(p) 



(5.3) 



for i — s,p (since there is no u structure in vacuum). The contributions to the microscopic side of the sum rule are 
given in the previous section. With the subsc ript labelling the contributions due to each figure, we find in the nuclear 
rest frame using the constraint in Equation (4.13): 
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3 2 7T 3 



,2 3 



2tt 3J 



M B E 2 L-* 



3% 5ys °M 2 

bn o 4 2a 2 , 



C j (Ml) = 0. 
Here we have used for the condensates 



a Q m 



1 



3 • 2 2 tt 7 M 2 



(5.4a) 



(5.4b) 



(5.4c) 



a 
a N 

bo 
b N 



- (2tt) z (qq) = 0.55 GeV a 
-(2tt) 2 (N\qq\N), 
(g 2 G ■ G) = 0.474 GeV 4 , 
(N\g 2 G -G\N) , and 



(qg s a ■ Gq) 



0.8 GeV 2 . 



The left-hand side of the sum rule is then 



We have taken into account perturbative 



corrections using the renormalization group anomalous dimensions of the interpolating field and the local operators 
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used in constructing the LHS. Each term is multiplied by the appropriate power of the anomalous dimension factor 



L 



= In (Mg/Aqcp) 

m (/V A QCD) 



(5.5) 



Since it arises from a conserved current, the anomalous dimension of the pseudo-Goldstone boson is zero. We set the 
renormalization point of the OPE p = 500 MeV and the QCD scale parameter Aqcd = 100 MeV. In the continuum 
model, the contributions from higher-energy states are included by using the prescription whereby terms proportional 
to M B n are multiplied by factors of E n _\ defined as 



" 1 

E n = l- exp {-ay M%) ]T j 



k=0 



k\ \ M 



We introduce a pair of functions $W as 



$W {Mis) =^>(Ml,.. 



A 1 



$W (M 2 ,s) 



2m 2 
Mj 



-mtt*) (Ml,s), 
o ( 2m 2 ~ 



e m b Am + me m b AA 2 — 2asoe m ~b A 



so 



M 2 B 



A 2 e m b Am - me 



-IT* 



Ai3 • 2 5 5 



^aa 2 

h 



Ai3 • 2 5 



2 T s 

a L9 s 
A13 • 2 2 



A13 • 2 2 



e M B As 



(5.6) 



(5.7) 



(5.8) 



The sum of these functions has no dependence on the residue shift AA 2 , and sufficient cancellation occurs between 
the terms proportional to the threshold shift Aso that we can ignore these terms. The in-nucleus mass shift of the 
nucleon can then be determined using 



1 



c 2 A 2 



Am = exp ( ^ ) *<•> {M% , s) + *« (Af 2 , S ) 



(5.9) 



We adopt a derivative improvement to the sum rules suggested by Drukarev and Ryskin pj| where 



to 2 dM 2 B 



(5.10) 



This improvement minimizes the effect of the factorization approximation of the four quark condensates by eliminating 
the terms proportional to p~ 2 from the dominant cont ribu ting diagrams, therefore aiding the convergence of the sum 
rules. Applying the derivative operator to Equations (|5.4|) we find 
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Re [I {a, c)] 
Re[/(a,c)] 
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(5.11c) 



The improvement does not remove all dependence of the LHS on the four-quark condensate, but these remaining four- 
quark terms are largely suppressed in comparison to terms with similar Borel weighings. As with the unimproved 
sum rules, a linear combination of these expressions can be built that has no dependence on the residue shift and 
negligible dependence on the threshold shift, giving the mass shift of the nucleon as 



Am 



1 



c 2 A 2 



— exp 



m 



4 s) (M%,s) -rm/> (p) (M%,a) 



(5.12) 



The improved sum rule is used in the analysis of the nucleon mass shifts in this work. 



B. Beyond the Chiral Limit 

As discussed above, in Ref. || it was shown that processes illustrated in FigureB produce terms proportional to 
m^// 2 . Also the m 2 Zn(m 2 ) terms were shown to cancel using the analysis of Ref. ]15[ |. It was thereby shown that our 
pion cloud model is consistent with chiral perturbation theory to leading order in the nonanalytic contributions to 
the nucleon mass in vacuum, and that we could use the chiral perturbation theory result to constrain the parameters 
of the model 




FIG. 6. The diagrams corresponding to the lowest-dimension chiral corrections to the nucleon correlator in medium. 



In dense nuclear matter (and at finite temperature, Refs 
seems to occur in the QCD sum rule method M. Equation ([1101 



29]) another problem with chiral symmetry breaking 
can be rewritten as 



= -p 



(5.13) 



However, the LNAC to the ttN sigma commutator is proportional to m\ |23|). The leading nonanalytic contribution 
to the in- medium quark condensate is thereby seen to be of order pm T , which is forbidden by chiral counting rules 
p?j| . Birse and Krippa accounted for this by adding contributions to the RHS where the nucleon current interacts 
with soft pions from the pion cloud of another nucleon in the nucleus Pq| . Using the soft pion theorem of PC AC 
pl| they find a term proportional to pm^ on the phenomenological side which cancels the offending term on the 
microscopic side arising from the in-medium quark condensate. As we go away from the chiral limit in our pion cloud 
model, we adopt the same modification to the phenomenological nucleon correlator, and in doing so we need to make 
sure that there are no contributions to the OPE from the pion cloud that would violate the chiral counting rules to 
leading nonanalytic order. This is indeed the case, since the chiral expansion of the pion propagator contributes lowest 
order terms proportional to m 2 , maintaining the leading order behavior. Likewise, similar to the vacuum contribution 
illustrated in Figure ||a lowest dimension in-medium contributions shown in Figures ||a and||b are possible. However, 
their effects are on higher order nonanalytic terms as they provide corrections proportional to pm\ and pvn\ ln(— m 2 ). 
Thus, to first order in baryon density the pion cloud model for the nucleon correlator in a nucleus is consistent with 
chiral perturbation theory in the leading nonanalytic contributions to the nucleon mass. 
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VI. RESULTS FOR THE NUCLEON MASS SHIFT IN SPIN-ISOSPIN SYMMETRIC NUCLEI 



This model is best applied to light nuclei where there are large surface volume regions and therefore the density 
gradients extend over a large fraction of the nucleus. In infinite nuclear matter these denisty gradients are zero, so 
that the s-wave part of the nA optical potential, which we neglect dominates the in-medium pion cloud contribution 
to the nuclcon correlator. In matter with Z = N, the strength of this potential is very small and the contributions 
of the diagrams in Figure || to the nucleon mass are negligible. As a result, the pion cloud has almost no effect on 
the nucleon mass in infinite nuclear matter, similar to the pion cloud effect on the nucleon mass in vacuum. It is 
evident that in large nuclei the pion cloud contribution to the nucleon mass becomes non-negligible only near the 
nuclear surface, a well-known property of the Kisslinger potential. So long as the nucleon remains a distance away 
from the surface comparable to about twice the ttN interaction range (with r^N ^0.7 fm) our mass shift result for 
this nucleon reproduce the work of previous authors Jl|,0 . Also application to large nuclei would require an extension 
of the method to account for N ^ Z. 

In this paper, we work with models of two small symmetric nuclei that arc familiar from ixA scattering e xperi ments, 
16 O and 40 Ca. The nucleon density distribution of the ground states of both nuclei are given by Eq (fufl ). For 
16 we choose the parameters c = 1.75 fm and a — 2.0 |32|]. This produces a distribtuion of nucleons with a local 
minimum at r = and an RMS radius of ~ 2.6 fm. The parameters for the 40 Ca density distribution are taken as 
c = 2.80 fm and a = 1.1. 
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FIG. 7. The mass shifts in the 16 nucleus at 0.5po with constant pion cloud probabilities. The solid line is the total mass 
shift, the dashed line is the mass shift due to the correlator without the pion cloud, and the dot-dashed line is the shift due to 
the pion cloud contribution. 
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FIG. 8. The mass shifts in the nucleus at po with constant pion cloud probabilities. The solid line is the total mass shift, 
the dashed line is the mass shift due to the correlator without the pion cloud, and the dot-dashed line is the shift due to the 
pion cloud contribution. 



We find the mass shifts in the chiral limit, ignoring in our analysis the effects of the small l eadin g nonanalytic terms 
discussed above. The mass shift of a nucleon in each nucleus is found by evaluating Eq. (5.12) in a Borcl window 
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0.8 GeV < Mg < 1.4 GeV. If the plot of the mass shift as a function of the Borel mass is relatively flat in this window, 
then the sum rules are sufficiently convergent and the mass shift is taken at this flat value. To test the stability of 
the sum rules we look at the mass shift at two points in the nucleus: first in the middle of the surface where the local 
density is half the nuclear matter density pq, then at the center of the nucleus, where the central density is equivalent 
to pq. As can be seen in Figures 0]l^, the sum rules in both nuclei provide sufficiently stable solutions in the Borcl 
window. 
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FIG. 9. The mass shifts in the 40 Ca nucleus at 0.5po with constant pion cloud probabilities. The solid line is the total mass 
shift, the dashed line is the mass shift due to the correlator without the pion cloud, and the dot-dashed line is the shift due to 
the pion cloud contribution. 
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FIG. 10. The mass shifts in the Ca nucleus at po with constant pion cloud probabilities. The solid line is the total mass 
shift, the dashed line is the mass shift due to the correlator without the pion cloud, and the dot-dashed line is the shift due to 
the pion cloud contribution. 
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FIG. 11. The mass shift from the correlator structures in the 40 Ca nucleus at po with constant pion cloud probabilities. 
The solid line is the mass shift from vector contributions to the correlator, and the dashed line is the mass shift due to scalar 
contributions to the correlator. 

The convergence is amplified in the plots of the individual contributions of the vector and scalar correlators (Fig. 
[TT| ). We find that the nucleon mass shift in an 16 nucleus is negligible at the central baryon density, and —150 MeV 
at one-half central density. The mass shifts in a 40 Ca nucleus are —260 MeV at p = po and —400 MeV at p = po/2. 

Some pion cloud effects might be related to the spontaneous breaking of chiral symmetry. Since the reduction in the 
value the quark condensate is a signal of partial restoration of chiral symmtery in the nucleus, the pion cloud amplitude 
C2 might decrease accordingly. We explore this idea by modifying our model, giving the pion cloud probability c\ a 
linear dependence on density such that it has the vacuum value at zero density, and it vanishes at the density where 
the quark condensate (qq) p goes to zero. Similarly, the probability of no pion cloud c\ increases from the vacuum 
value of 0.5 to the phase transition point value of 1.0. Figure fl2| shows the mass shift of a nucleon in a 40 Ca nucleus at 
central density utilizing this modification to the model. The shift in 16 O shows a similar, if less dramatic, reduction. 
As in the constant pion cloud probability calculation, the sum rule calculation is sufficiently flat in the Borel window 
to provide a value for the mass shift. The mass shifts in 16 are ~ MeV at central and —100 MeV at surface 
density. In 40 Ca we find -100 MeV and -275 MeV, respectively. 
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FIG. 12. The mass shifts in the 40 Ca nucleus at po with pion cloud probabilities linearly scaling with density. The solid line 
is the total mass shift, the dashed line is the mass shift due to the correlator without the pion cloud, and the dot-dashed line 
is the shift due to the pion cloud contribution. 

We compare these two methods to the sum rules without an explicit pion cloud throughout the nucleus in Figures 
|l3| and [l4|. In these calculatio ns, t he density used at each point in the sum rules is the local density obtained from 



the density distribution (Eq. (3.16)) and the nucleon mass is plotted as a function of the distance from the center of 



the nucleus. Each mass shift is evaluated at the Borel mass Mb = 1.0 GeV. 
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FIG. 13. The mass shifts in the 16 nucleus as functions of distance from the nuclear center, calculated at a Borel mass 
of 1 GeV. The dashed line is the density distribution. The solid line is the nucleon mass with linearly scaling pion cloud 
probabilities, the dot-dashed line is the mass with constant probabilities, and the dotted line is the nucleon mass without the 
pion cloud contributions. 



In an 16 nucleus, the mass shift due to the pion cloud is small compared to the mass shift from the non-pionic 
terms but positive out to a distance of about 4.5 fm from the nuclear center. As the nucleon appraoches the nuclear 
surface the pion cloud contribution to the mass shift becomes positive and larger than the mass shift without pion 
cloud. Even though the local density in the interval 3.0 fm < r < 4.0 fm is small (hence the small mass shift due 
to the correlator contributions without the pion cloud), the pion cloud is still sampling a significant range of density 
gradients leading to a sizable optical potential strength, and therefore a large mass shift. The two different approaches 
to the pion cloud amplitudes in 16 O do not yield significantly different results throughout the nucleus. In the 40 Ca 
nucleus, the mass shift due to the pion cloud is always negative and comparable to the non-pionic mass shift. As with 
the oxygen nucleus, the largest mass shift from the pionic terms occurs near the surface, approaching —600 MeV. 
The constant and linearly decreasing pion cloud probability calculations converge to the same value in this region. 
The nuclear interior shows a significant difference between the three methods. If ci and ci remain constant over the 
range of densities, then the pion cloud provides a rather large mass shift (~ —275 MeV at central density). Once Ci 
and C2 are dependent on density, the pion cloud contribution is more than halved at the center of the nucleus, with a 
mass shift closer to —100 MeV. Clearly the choice of pion cloud method has a significant effect in larger nuclei. 




FIG. 14. The mass shifts in the Ca nucleus as functions of distance from the nuclear center, calculated at a Borel mass 
of 1 GeV. The dashed line is the density distribution. The solid line is the nucleon mass with linearly scaling pion cloud 
probabilities, the dot-dashed line is the mass with constant probabilities, and the dotted line is the nucleon mass without the 
pion cloud contributions. 



VII. CONCLUSIONS 



We have described a model based on QCD for the nucleon in a finite spin and isospin symmetric nucleus including 
a pseudo-Goldstone boson component of the correlator This approach is based on the model proposed by one of us in 
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Ref. H , with the extension of the pion free propagator to the in-medium propagator using an optical potential which 
is first order in the nuclear density. A fit to chiral perturbation theory restricts the parameters of the pion cloud 
component. Our main conclusion is that there are significant contributions to the mass of the nucleon in the nucleus 
arising from the pion cloud component of the nucleon. 

In the deep interior of the oxygen nucleus we find a small positive mass shift which becomes large and negative 
near the surface. In the calcium nucleus, however, the mass shift is everywhere large and negative. We introduce a 
parametrization of the pion cloud probability c\ which decreases linearly from the vacuum value to zero at the chiral 
symmetry restoration as a possible test of partial restoration of chiral symmetry. At high densities, the effect of the 
pion cloud is then reduced leading to smaller mass shifts. 

One obvious extension to this model is the consideration of the strange sector. By allowing for the inclusion 
of kaon cloud terms in the baryon current, we can construct similar sum rules for the mass shift of hyperons. The 
primary obstacles are the uncertainties of several in-medium matrix elements on the microscopic side, the complicated 
structure of the dispersion relation on the phenomenological side |24|], and the difficulty in finding an appropriate 
form for the K-A optical potential to determine the kaon self-enegry. The model can also be modified to handle N 
7^ Z nuclei. The nucleon correlator in such nuclei without the pion cloud has been discussed in Ref. and the 
optical potential is easily treat ed in asymmetric nuclei by including the isospin-dependent and spin flip terms in the 



nN scattering amplitude (Eq. (3.8)). A third venue of further research concerns the applications of the model to sea 
quark distributions in nuclei and the development of an improvement to QCD calculations of nuclear deep inelastic 
scattering in the interval 0.1 < x < 0.3. 

We would like to thank Mountaga Aw and Otto Linsuain for helpful discussions. 
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